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BACKGROUND 

Pulse outers typically measure and display 
various blood flow characteristics including but not lL* „ 
t. Moo* oxygen saturation of hemoglobin iTartriat lilT 
volume of individual blood pulsations and th^elf b^ 
Rations corresponding to each heartbeat^* LTpItient 
The oximeters pass lioht tto-^h _ Patient. 



a finger, an ear, the 
illy sense the 
The amount off 



«u , • • "» ww nearajeat of the nation*- 

*he oximeters pass light through human or animal bc^ * 
where blood perfuses the tissue such — - - 
nasal septum or the scalp, and photo. 

change in absorption of light in the tissue. The - 
light absorbed is then used to calculate^ amoirT^ 
constituent being measured. ° f blood 

be 0 f m WBTOd the tlBS »* i« selected to 

be of one or »o.re vavelengths that is absorbed by the Wood in 
an amount representative of the amount of the blooT 
constituent present in the blood. The amount of teansmi^ 

»SrC[r °~ nt * - tissue J^a 

reaches the light detector. Another sourc7lTL^ 
alectroaagnetip coupling from ctheTel^! T T ^ 

r r- Hotion 2 «• - * ^trr^ 

For Katnc, «o»Uig, thQ coupllno h "l ™ , ' 
£t«>o. In .« Itl o», .inc. „ M U „ Ilala , n ^ l f0t 
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resulting in a momentary change in volume at the point the 
oximeter probe is attached. 

I****.! - aotlon °** degrade 8i9nal "»* *« 

medical decisions, with the clinician being unaware of it. 

ZlJZ T Cl ^ y T " 18 rea0te storing of the 

patient, the motion is too email to be observed, the clinician 
is watching the instrument or other parte of the patient and 
not the sensor site, or in a fetus where motion is hiMen 

In one oximeter system described in U.S. Patent No 
5,025,79!, an aocelerometer is used to detect motion. «he„ 
motion is detected, readings influenced by motion are either 
eliminated or indicated as being corrupted, in a typicaT 

, meaSUreaontS ***** ^ the peaks and valleys of the 
blood pulse signal are used to calculate the desired 
characteristic. Motion can cause a false peak, resulting in a 
measurement having an inaccurate value and one which is 
recorded at the wrong time, in 0.8. Patent Ho. 4,802,486 
assigned to Helicor, the disclosure of which is incorporated 
herein by reference, an EKO signal is monitored and correlated 
to the oximeter reading to provide synchronization to limit 
the affect of noise and motion artifact pulses on the oximeter 
readings, this reduoes the chances of the oximeter locking on 
to a periodic motion signal, still other systems, such as 
that set forth in U.S. Patent Ho. 5,078,136, assigned to 
Helicor, the disclosure of which is incorporated herein by' 
reference, use signal processing in an attempt to limit the 
affect of noise and motion artifact, the .«6 patent, for 
instance, uses linear interpolation and rafce-of change 
techniques to analyse the oximeter signal. 

-i*-, *** ntttUra ° f readings impose a number of 

difficulties in dealing with noise. The oximeter re^eTon 
mathematical analysis of the reading, at two dtff erenu 
wavelengths. Because different amounts of light are absorbed 
at each wavelength, the magnitude of the motion artifact duT 
to the same motion will be different for each signal. This" ia 
complicated by the fact that the lights are alternately 



CONFIDENTIAL 
I ATTORNEYS EYES ONLY 



M 2412 

JA 29090 



Copied from 10 7 79033 on 07/1 5/200$ 



WO 97/00041 



PCT7US96/10296 



pulsed, and thus each is influenced by a different amount of 
notion, since the notion varies with tine. 

One system, set forth in PCT Publication Ho. 
WO 92/15955 (Vital Signals, inc.) correlates the non-noise 
portion of two wavelength signals and generates a noise 
reference signal. The noise reference signal is then provided 
to an adaptive noise canoeler to eliminate the noise from the 
desired signal. • 

Patent Ho. 4,714,341 discloses the use of three 
different wavelengths, rather than two, in order to detect 
when nbise is present. This patent teaches using the first 
and second wavelength signals to produce a first oxygen 
saturation value, and then using the first and third 
wavelength signals to produce a second oxygen saturation 
value. The two calculated values are then compared, if the 
values are equal, as they should be absent motion, the signal 
is presumed to be good, if the valuee are different, the 
signal is assumed to contain motion and is disregarded, 

SUMMARY OP THB INVENTION 
The present invention is based on analysis of the 
signal intensity received by the detectors, without separately 
measuring the motion signal, without providing feedback to 
cancel the motion signal and without attempting to 
mathematically eliminate the motion signal. Instead, the 
present invention mathematically recognises the presence of 
the motion signal and recognises a few key characteristics of 
the motion signal. First, although the magnitude of the 
effect of motion on the signal intensity for each wavelength 
will be different, the change in the logarithm of the motion 
component will be approximately the same (for signals at 
approximately the same time) . This allows the motion 
component to be cancelled out in a ratiometrio equation. 
Second, it is assumed that the blood pulse signal is not 
affected by motion. This second assumption is more of an 
approximation, since the blood pulse signal is somewhat 



I MA8105162 

.^CONFIDENTIAL 
! ATTORNEYS EYES ONLY 



M2413 

JA 29091 



Copied from 10 7 79033 on 07 / 1.5/ 200 S 



WO97/D004I 

™* notion*. ^ „, tUfl'T *-* not 
don. ,or Ittor . ^ wvuI^TL" S!" ™" °" *• 

.awlMgth =130.1.. The lo,«lth» „ ^L"!! f "*™ 0a 

» TO ""^"tzTt! .TIT " tlM iB 

-^."on «Ooh i. . roootio.™ rj^^* "™« » :- 

— o. s :ir;.r^Lr„:rr i? 

assumed to ha «.>,„ .... . - 18 with tltte ana is 



f MAS 100163 

!atto < B^ only 
M2414 



JA 29092 



Copied from 10779033 on 07/1 5/200$ 



WOW/00041 



PCT/US9M0J96 



•unfi* stalled description taken to conjunction with the 
accompanying drawings. • 

BRIEF DESCRIPTION OP TOE DRAWINGS 
- Pigs, ia-id are diagrams, of an intensity signal 
showing the effects of puisatile flow and motion noise; 

Pigs. 2A and 2B krc diagrams illustrating the effect 
of motion on the path length of emitted light, and thus on S 
intensity of received light; and 

Pig. 3 is a block diagram of a system according to 
the present invention. 9 

DETAILED DESCRIPTION OF THE' PREFERRED EMBODIMENTS 
Pigs. 1A-1D illustrate aspects of a pulse oximeter 
signal which the present invention. takes advantage of 
Pig. 1A shows the logarithm of a detected infrared signal ' 
Pig. IB shows the logarithm of a detected red wavelength ' 
signal For both of these figures, the signal includes motion 
occurring in the interval of 5-12 seconds. Otherwise, both 
the red and infrared signals are noise-free optical signals. 
Pig. 1C shows the result of a subtraction between the signals 
in Pigs. IA and IB. As this illustrates, the subtraction 
cancels out the noise. This is because the data exists in 
logarithm form, and the motion corruption is additive. 
Accordingly, in addition to calculating saturation the 
difference waveform (Pig. ic) can be scaled, and then 
^btracted from either the logarithm of the IR or red signal 
«t£a£? ^ * 8tlaate ° f B0tlOn "> "*owTthis 

Pig. 2A illustrates one possible example of how 
action can effect the intensity signal. A light emitter 16 is 
shown emitting rays 18 through a patient's finger 20. This is 
detected by a detector 22. As can be seen, the distance tJ« 
the emitter to the detector, 0, will determine th.^ount e 7 
light emitted by the emitter reaching the detector, since . 
there will be a natural spreading affect of non-collimated 



MAS 106164 
, m CONFIDENTIAL 
I ATTORNEYS EYESONLY 



M 2415 

JA 29093 



Copied from 10 7 79033 on 67/15/200$ 



WO97/D0WI 



PCT/US96/10J96 



light rays. The farther away the detector is, the nore 
spreading results. \ 

Fig. 2B illustrates another example showing how 
motion of a finger can compress and widen the finger 
(exaggerated in the figure) and temporarily cause the light 
emitter X« to move away from the detector an additional 
distance indicated by arrow 24. This additional distance will 
cause less of the light to reach the detector, since there 
will be more spreading of the light emitted at this larger 
distance. This will result in a lower intensity waveform 
being detected by the detector. Alternately, compression 
could result in a higher intensity waveform. Motion and noise 
can take other forms as well; and can vary for other reasons 
than non-ccllimated light rays. For instance, the emitter and 
detector could be slightly misaligned. ' 

The present invention recognizes that the 
calculation for determining oxygen saturation by pulse 
oximetry using .the -ratio of ratios- can be assumed to have a 
motion term Which is independent of any particular wavelength 
An understanding of this first requires an understanding of 
how the ratio of ratios is calculated. 

Using Laabert-Beeris law as a starting point, 
equation (1) below is used to determine saturation in pulse 
oximetry: 

XU, t) - X,U>exp<-<sB 0 U) + (x . fl , ff(X)) J(t)) <*> 
where t X a wavelength 

t " time 

X 0 - intensity of light transmitted ~ — * 

X - intensity of light detected 
a o oxygon saturation 
Pp< P, - empirical^ derived absorption coefficients 
fox oxygenated end deoxygenated hemoglobin 
respectively ' 
It) • a combination of concentration and path 
length from emitter to detector as 
a function of time 



AT T°jWEYBEYE80NLY I 

M2416 



JA 29094 



Copied from 10 7 79033 on 07 / 1.5/ 200 S 



tr * dltional ^ to solve for ratio ot 

ratios and then calculate saturation. 

Take natural logarithm of equation (1, for IR ana Red, 

103Xel °9 X o-<Spc*(W)P r )i (2) 
Differentiate equation (2) yith respeot to time: 

Divide Red (3) by IR (3) 

dlogj (X.) /dt gPo (1 . g) fl _ ( , 



For a discrete time sample, equations of the above form M „ k. 
rewritten by noting: Y * f0rB 0311 ^ 



Using logA-logB - logA/B, the above equation c 
written as: 



Where r io the "ratio of j 
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solving (4) for a using (5) gives: ; « 

FroB (5) note R can be calculated using two points 
corresponding to measurements at two different tines, t. 
Alternately, a family of points can be used. 
To see this latter point define: 

Then, equation (5) can be written as: 
T(t) - xx(t) 

and for a family of points over time this will describe a 
cluster of points that define a best-fit line of y versus x 
with a slope given by R. 

The present invention modifies the above equations 
by recogniaing that a tern can be added to account for notion 
and noise, m particular, the motion and noise coaponent can 
be represented by a function which varies with time and is 
wavelength-independent. This recognition allows a 
mathematical solution to isolate and eliminate the motion and 
noise components without requiring prior art methods such as 
separately measuring the motion. , 

Hatiaa. For example, to account for motion and 
noise, we can modify equation (X).by multiplying by a time 
varying function « (t) representing wavelength-independent 
motion or noise. This gives the following equation: 

ju,t> - x»a)n(t)exp(.(sp 0 a) + (i- S )p r a, )J(t)) (6 , 
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We can then solve for s using the sa »e steps as used above. 
First, we take the logarithm! 

log J-logJ < , + log 1 ,.( a p jf+ a . ff)ppi 

s * 
Next, we differentiate with respect to time: 

^ - -«p -.<•».♦ c M w || .'; 

Then ' we determine the ratio of Red to XRt 

^ if d ipg , /dt is large ^ ^ other ^ 

ratio of ratios will be driven towards unity, driving 7 
towards a wavelength-dependant constant. So because In «h„ 
aodel optical coupling due to motion appears ide^UcaSyto 
both wavelengths, its presenoe drives the saturation^ taTs 
wavelength-dependant constant. 

*** PreSent ±nVention allows a calculation of 
blood oxygen saturation by mathematically recognising^ 
»otion signal. This enables a solution vhichdoesTot^L^ 
separately measuring the motion signal, proviCfe^cT* 
cancel the motion signal, or attempting to .nathemaiicafT 
eliminate the motion signal. Set forth below are *T 
preferred embodiments for implementina th« , 

A Thgflft-waval enert ,h ff nl^ t lff n 

Let ^ be some other wavelength (not IR or Red) 
How tajca the logarithm and differentiate this thirl 
wavelength, obtaining (7) . Gne approaoa ^ ^ J. 
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difference 1R with this new wavelength, and similarly with 

Here is a batter solution. Rewrite (7) as: 
How to introduce, some matrix algebra, define: 



*-^logij 
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With this notation * 



k» «fc Ha) 

b*Cx 

JB » X) 



a. " long „c 1. run r„o,, Uf 

saturation. 1 ** £or » 9iven 

c new wavelength is one that oivaa th« >,j„. 
*»0 option c=.«io*^. t „ th,^ t^lT^- 
— ««. In no notion, th. ..t^nt". oSon^^STt 1 "- 
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With two wavelengths we have: 



[CU °22 I, 



to Mtu egUaWons and «*nowns. one approach ls 

to return to calculating J? by rewriting (8) : 8 



• where m l s the notion tern, aa defined earlier 

* ratios, ana v is the signal with no motion!' "* 
There are two key assumptions which make the 

?*l7o?:*2T*' tixsti altho °* "» ofL 

STeLT < ~ °° lntenai ^ 8i9nal * Cerent, 

the change in the logarithm of the motion component at twT 
d liferent times will be the sane (which assumes the Sr^ent 
time signal sample, are adjacent or close together in^T 
This allows the motion component to be cancelled^? £T ' 

l277o° T ti0ft ' 6eP0nd *• that Se 

motion does not cause any effect on the remainder of the 
Ration. There is sone effect, since motion Z change the 
pulse flow characteristics of the blood, but this is"»ic^w 
a small effect compared to the motion when there U °^ 
n^f 7 ° nPre8Bnt - that the motion has 

T ■ ^ 0t °^*»*ion measurement^ 

assume that v and a are not related. 

Another assumption is that the amount of motion la 
the same at the time of both intensity signal measured ! 
^wavelengths. This is a iJL^^SE^* . 
^-typical motion signal varies at., rate of arounTa £ 
* a he light pulsing frequency is typically atTrate of 
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Assuming v and m are independent over time, ty ffl ,» 0 
fO^scae defined inner product. Substituting fer y and m 

Solving for the R that solves this equality yields. 



There are two problems with this approach. 

Sbov* ^ * aPPr ° aChe3 ° ne ' ^roaches aero, and the 

above equation approaches zero divided by aero. This f a m- «? 

z ^.r ~ tor *"•" ^ - ~ - 

A acre limiting problem is the assumption that 
fv,»,=o certainly the motion signal is independent of the 
arterial pulsatile signal, but during motion, v also has^ath- 
length concentration effects in it that are highly correlate 
with », thus biasing * away from its true valuT 

Pig. 3 is a block diagram of one embodiment of a 
^oximeter implementing the present invention. Light from 
t*Ds X4 passes into patient tissue «, and after being 
transmitted through or reflected from tissue «, the light is 
received by photosensor «. Either two or three iW can be 
»aed depending upon the embodiment of the present invention! 
Photosensor 16 converts the received energy into an e^trLl ~ 
signal, which is then fed to input amplifier 20. 

tight sources other than LEDs can be used. For 
ute7ir: h la8erS T" U8e *' ^ ° "™» «uld be 

r: prlata mters eith ~ - - 

Time Processing unit (TPO) 48 sends control 
signals 68 to the led drive 32, to alternate* activate the 
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LEDs * ****** depending on the embodiment, the' drive nay 
control two or three LED s. \. 

«he signal received from input amplifier 20 is 
passed through three different channels ae shown in the 
embodiment of Pig. 3, for three different wavelengths. 
Alternately, two channels for two wavelengths could be used. 
Bach channel includes an analog switch 40, a low pass 
filter. 42, and an analog to digital (A/D) converter 38. 
Control lt**s 69 from **u 48 select the appropriate channel at 
the time the corresponding LED 14 is being driven in 
synchroniMtion. A queued serial nodule < 0 SH) 46'receives the 
digital data from each of the channels. CPU so transfers the 
data from QSH 46 into HAH 32 as Q8K 46 periodically fills up. 
to one embodiment, QSH 46, TO* 48, CPO SO and BAM S2 are part 
ot one integrated circuit, such as a DMC68HC16 microcontroller 
from Motorola. 

She method of the present invention is practiced by 
CPU 50 on the data in RAM 52 as received through the various 
channels from photodeteotor 16. -me signal from 
Photodetector.16 is the signal which originated from LEDs a4 
as reflected or transmitted by patient 18, and including • ' 
undesired noise artifact. 

As will be understood by those of skill in the art 
tha present invention can be embodied in other specific forms 
without departing from the spirit or. essential characteristics 
thereof. For example, saturation could be determined using 
different mathematical calculations, once it is recognised 
that the motion term is a function of time that is independent 
of wavelength and is approximately the same for two adjacenT 
in time signal samples at two different wavelengths, m one 
example, the mathematical determination could be done by 
dividing the two intensity, equations to eliminate the motion 
term. Although this would only eliminate the motion from one 
wavelength equation, this could be done for alternate 
wavelengths in alternate samples. In a three wavelength " 
embodiment, division of two separate pairs could be done to 
eliminate the motion signal. Accordingly, the disclosure of 
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Siu^f emb ° diffiBnt ° f intent tQ ba 

^< ? «- -pa of the inventL 
wnion j.8 set forth In the following claims. 



' ATTORNEYS EYE8 ONLY. 

M2425 

M 29103 - 



Copied from 10 7 79033 on 07/1 5/2005 



WO 97/00041 



PCT/DS9M0396 



wavelengths; rent 

determining said saturation by manipulating 
Said ^ - d density signals^^ * 

assumptions that ; 

i) an amount of notion is the same at *h- 
mm ti_ ,„ eaeb of „u i^ity .i^/LT 

ii) the notion components of said 
intensity signals are proportional to one another. 

2. The method of claim 1 wherein said , 
*tep assumes that the derivative of the^ ^TSe ^ 
motion components o* said intensity signals are the same"! 

steps of, 3 * Beth0d ° f Clala 1 fUrth - -"Prising the 
representing each of said intensity B i ma1 
• function of said saturation, the ' " 

corresponding to the intensity signal, and a time- 
variable motion term corresponding to motion 
motion terms being proportional tl ' 
or said intensity signals, ^ '° r ° ach 

taking the logarithm of each rebr MM *-.« 
said first and second intensity ai^l.? ** °* 

differentiating each logarithm, 
equating the first differentiated 
the first intensity ,ig„ al to v + m, SST.^TST * 
Portion pf the 8ignal toe to aotion ; «■ » is «"» 
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equating the second differential:."! ««. 
the second intensity signal to Rv TT^Zt IV * 
of first and second wavelength ratios^ I lZ 
ratio being the logarithm of the ratio „ ^ 
signal for the length at I£ 

expressing said representations as a matrL, 
solving said matrix for r h V 4 

determining said saturation from R. 



4. 



prodatonained _ t f „ . p^^.t.^!'*" * 

wavelengths; "terent 
„„„ ■ sensin 9 **• intensity of said radiation for 

Wil T T elen9thS ^ At P — througTa 
portion of said patient to produce first «~,Z. . 
third intensity signals; ' 8eC0nd "* 

•representing each of eaid intensitv «i m „ 
afnnction of-said saturation, the vave^S *° 
corresponding to the intensity signal, andVtime- 
variable motion term correepondCto mo^on no! e 8ald 
»otion terms being proportional to one anoth^r'eact 
of p«id intensity signals; and * 
solving the three functions to obtain a value 
for said Saturation. ™ lue 
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7. The method of claim i wherein said deterging 
step assumes that the derivative of the logarithm of the 
motion component of each intensity signal is the same. 

8 . the »ethod of claim « wherein each of said 
functions includes a plurality of coefficients, and further 
ocmprising the step of determining a set of ooef f icients f or 
said third intensity- signal from a measurement in the VJZe 
of motion noise and a determination of said saturation froT 
said first and second intensity signals. ' r °* 

9. The method of claim 6 further comprising the 

stops oft 

eal - « . taki " 9 tl18 logarltto of "Presentation of 
said first, second and third intensity signals; 

differentiating each logarithm; 

putting the differentiated logarithms into a 
matrix; and 

solving said matrix for Baid saturation. 

10. & method for measuring the saturation of a 
blood constituent in a patient comprising the steps of j 

irradiating said patient with electromagnetic 
radiation of three discrete, different wavelengths; 

sensing the intensity of said radiation for 
each of said wavelengths after it passes through a 
portion of said patient to produce first, second and 
third intensity signals; 

representing each of said intensity signals as 
e function of said saturation, the wavelength 
corresponding to the intensity Signal,' and a time- 
variable motion term corresponding to motion noise, said 
notion term being the same for each cf said Intensity 
signals; and * 

solving the three functions to obtain a value 
for said saturation. 
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1 ' lx - »ethod,of claim X0 wherein each of said 

2 functions includes a plurality of coefficients;- and further 

3 comprising the step of determining a set of coefficients for 

4 said third intensity signal from a measurement in the absence 

5 of motion noise and a determination of said saturation from 

6 said first and second intensity signals. 

1 12. The method of claim 10 further comprising the 

2 steps oft 

3 taxing the logarithm of each representation of 

4 said first, second and third intensity signals, 

5 differentiating each logarithm; 

6 putting the differentiated logarithms into a 

7 matrix; and 

8 solving said matrix for said saturation, 

1 13. A method for measuring the saturation of a 

. 2 blood constituent in ff patient comprising the steps of: 

3 irradiating said patient with eleotrottagnetio 

4 radiation of two discrete^ different wavelengths, . 

5 sensing the intensity of said radiation for' 

6 each of said wavelengths after it passes through a 

7 portion of said patient separately to produce first and 

8 second intensity signals; 

9 representing each of said intensity signals as 

10 a function of said saturation, the wavelength 

11 corresponding to the intensity signal, and a time- • 
J? ^f 19 ttotlon *■» corresponding to motion noise, said 
is motion term being the same for each of said intensitv 

14 signals; 1 

" taW »* toQ lo**ithm of each representation of 

16 said first and second intensity signals; 

17 differentiating each logarithm, 

" eguating the first differentiated logarithm of 

» the first intensity signal tov + l( » ls ^ . 

20 portion of the signal due to motion; • 
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- second ™ W> 

aetermining eaid saturation from R, 

J w 14 # *» W^atus for measuring the seta—** 

* blood constituent in a patient uprising, 8atu * atio ° <* a 

4 " flrSt 8ee °na emitters, said emitters 
J. emitting radiation of first ana second di«erenT 

5 wavelengths; •"went 

7 a detector >* the intensity of said 

U f V 8aid Mounted relative 2 !aid1L. 

it passes through abortion of said patient; 

a controller for alternately aotivatina 
emitters so that said detector detects the dirSent 
vavelengtha at diTferent times to produce fS"^ 
second intensity signals, and an<1 

control means for determining eaid saturat* „„ 
with the assumptions that y ^ La 

i) an amount of motion is the same a t «,« 

same time for each of said intensity signals, a^r 

. ii) the motion components of said 
intensity signals are proportional to one another. 

15. The apparatus of claim 14 wherain » «'« ^ 
means further comprises means for ImJ^^?^ 7*** 
of the logarithm of the motion oompon^^eW 
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J 6 ' apparatus ef **** " herein said 

jztz a puiss «* ~ means ^ 

Blona1ra aeaas for Weeentiag each of said intensity 
signals as fa function of Bal d saturation, the wavewL 

ZEST. t0 lntSnaity 8l *»l< *»* a t™ 91 * 

variable notion ten. corresponding to notion noise aaia 
mooter* being the same for each of seid in^nsV* 

means for taking the logarithm of each 
si™" 1 " ^ **' B6COnd **— * 

means for differentiating each logarithm 

means for equating the first differentiated 
logarithm of the first intensity signal to v Tl lt\ 
ia the portion of the signal due to Motion, ™ 

means for equating the second differentiated 
logarithm of the second intensity signal to ZTt ^ 
» ia a rati, of first and second ratios" ^7 

wavelength ratio being the logarithm of £e ratio of ^ 
intensity signal for the wavelength at first 1 set£ . 

matrix; ** Bii * "presentations as a 

v and » f0r /° 1Vln? 8aid ***** *>r R by assuming 

vend » are independent for some defined inner product 

jaeans for determining said saturation from *. 

17. The apparatus of claim 14 further * . 
diaplay, coupled to said control neans t^lZr 7 ^ * 
saturation. ' dls I**yi«g said 

18. The apparatus of claim 17 *,~*.v. 

alam coupled to said control .sensor iSS^T^" 
saturation is less than a predetermined emouT£r 7 
predetermined period of time. 
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«»e apparatus of elai» x 4 further comprising 

with eieetromagnetioradiation of a third discrete 
different wavelength; ' 

8aid COntroUer Alternately activating said 
first, second and third emitters so that said detector 
Produces first, second and third intensi^ signets 
aaid control means including 
means for representing each of said intensity 
signals as a function of said saturation, the 
wavelength corresponding to the intensity signal, 
and a time-variable motion term eorrespondinTto 
action noise, said motion term being the same for 
each of said intensity signals; and 

^means for solving the three functions to 
Obtain a value for said saturation. 
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